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ABSTRACT 


Selecting  the  best  Gamma  population  from  a  given  set  of  Gamma  populations 
is  treated  from  a  decision  theoretic  point  of  view.  Cost  of  sampling  and  pen¬ 
alties  for  wrong  decision  play  a  role  in  the  determination  of  optimum  common 
sample  sizes.  Minimax  sample  sizes  are  determined  under  two  different  penalty 


functions 


•  4 


’“i 

4 


& 

M 


AMS  1980  subject  classification:  primary:  62F07,  secondary:  62C20. 

Key  words  and  phrases.  Selection  problems.  Gamma  distribution,  cost  of  sampling, 

penalty  functions,  minimax  criterion,  minimax  sample 
size. 


$ 

,s*u 


ON  THE  SELECTION  OF  BEST  GAMMA  POPULATION 
DETEBMINATION  OF  MINIMAX  SAMPLE  SIZES 


r«rv  ✓ 


N.  Dailami 
M.  Bhaskara  Rao 


Department  of  Probability  &  Statistics 
The  University 
Sheffield  S3  7RH 
UK 


and 

K.  Subram anyam 


Department  of  Mathematics  &  Statistics 
University  of  Pittsburgh 
Pittsburgh,  PA  15260 


1.  Introduction 


Let  ir^,  7r2”***TTk  be  k  Gamma  populations  with  unknown  scale  parameters 


01,02 . 0fc  respectively  and  a  common  known  shape  parameters  c  >  0.  The  probability 


density  function  of  the  i-th  population  it^  (1*1,2,...  ,k)  is  given  by 


fjXx)  -  9~C(T  Cc))"1  exp  C-  *jp>xC_1,  0  <  x  <  <*>. 


The  main  objective  of  this  article  is  to  describe  a  method  of  selecting  that  Gamma 
population  with  the  least  0-value.  This  problem  has  considerable  bearing  on  the 
problem  of  selecting  the  best  of  several  normal  populations  in  the  sense  of 
selecting  that  normal  population  with  least  variance.  Basically,  any  selection 
problem  consists  of  two  components. 

1.  Draw  a  random  sample  of  size  n^  from  population  it^,  i*l,2,...,k. 

2.  Suggest  a  statistical  procedure  R  which  once  the  data  are  given  clearly 
spells  out  the  best  population. 


2 


In  many  selection  problems,  a  natural  statistical  procedure  R  manifests 
itself.  The  real  problem  is  the  selection  of  the  sample  sizes.  The  sample  sizes 
are  determined  following  certain  optimality  criterion.  Following  the  lead  given 
by  Somerville  (1954)  and  Ofosu  (1972) ,  we  determine  the  sample  sizes  taking  into 

-4 

account  the  cost  of  sampling  and  penalties  imposed  when  wrong  decisions  are  taken. 
The  nature  of  the  determination  of  the  sample  sizes  is  decision-theoretic  in 
character  and  we  adopt  the  minimax  criterion  as  the  optimality  criterion.  Ofosu 
(1972)  has  studied  this  problem  under  a  particular  penalty  function  and  determined 
the  optimal  sample  sizes  using  minimax  principle.  But  his  arguments  have  a  gap 
and  it  is  debatable  whether  the  sample  sizes  he  has  given  are  really  optimal 
sample  sizes. 

In  this  paper,  we  consider  two  types  of  penalty  functions  one  of  which  is 
slightly  less  general  than  the  one  introduced  by  Ofosu  (1972).  Under  these  two 
penalty  functions,  we  minimize  the  resultant  loss  functions  over  the  entire  para¬ 
meter  space  for  every  fixed  common  sample  size  n.  The  solution  to  this  minimiza¬ 
tion  problem  works  out  quite  explicitly  overcoming  the  gap  present  in  the  paper 
by  Ofosu  (1972). 

For  a  general  introduction  to  ranking  and  selection  problems,  the  reader 
may  refer  to  Gibbons,  Olkin  and  Sobel  (1977)  or  Gupta  and  Panchapakesan  (1979). 

2.  Preliminaries 

If  X^,X2,...,Xn  is  a  random  sample  of  size  n  from  a  Gamma  population  with 

n 

scale  parameter  3  and  shape  parameter  c,  then  £  X .is  a  sufficient  statistic  for 

i-1  n 

9  for  known  c.  Let  T  ■  T(X. ,X. ,... ,X  )  ■  n  T  X . .  Then  ET  ■  c9  and  T  has  a 

“12  n  i-i  i 

Gamma  distribution  with  scale  parameter  6n  and  shape  parameter  cn.  For  the  problem 
mentioned  in  the  introduction,  suppose  that  we  wish  to  draw  random  samples  of 
same  size  from  each  population.  Let  X^ ,  j  -l,2,...,n  be  a  random  sample  of  size  n 


from  population  i-l,2,...,k.  Let  ■  T(X  •  •  •  »X^n> »  i“l»2,...,k. 

Now,  a  natural  statistical  procedure  manifests  itself  for  selecting  the  best  Gamma 
population. 


Statistical  Procedure  R 


Declare  population  ir 


1.2. . . . .k)  to  be  the  best  i£  <  T  for  all  1  ^ i. 


Cost  of  Si 


We  assume  that  the  cost  C(n}  of  drawing  a  random  sample  of  size  n  from  all 
populations  to  be  of  the  following  form 


C(n)  -  c  +  c.n  ,  n«l,2,..., 

o  i 

where  cq,c^  and  d  are  non-negative  constants.  cq  and  are  measured  in  the 

same  units,  and  c  represents  fixed  administrative  costs  involved  in  setting 
o 

up  a  sampling  plan.  If  d  ■ 1,  the  cost  of  taking  additional  samples  rises  linearly 
with  n.  If  d  < 1,  the  rise  in  the  cost  does  not  increase  relatively  with  increasing 
sample  sizes.  If  d  > 1,  it  will  become  more  and  more  expensive  to  take  additional 
samples. 

Wrong  Decisions 

Let  9i»S2»* •• ke  a  configuration  of  the  parameters  of  the  k  populations 
^1*^2’ * ** ,7rk  resPectively*  Let  0q^  L 0(2}  —  •••  —  ®qc)  be  the  ordered  arrangement 


of  elfe2,...,ek.  l« 


, tt  correspond  to  the  parameters  0,n,9,»v 

(i/ 


respectively.  We  follow  the  usual  convention  that  if  two  or  more  0-values  coincide 
the  corresponding  i^'s  are  taken  to  be  in  increasing  order  of  magnitude.  Let 
T.  correspond  to  ir  for  j  -  ] ,2,...,k.  According  to  the  configuration 

J  j 

01,®2* *  *  * *®k*  wi  t*le  kest  population.  There  are  k-1  different  ways  of  going 
wrong.  These  are  listed  below. 
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Type  of 

wrong 

decision 


Description 


ir  is  declared  to  be  the  best. 
i2 


ir  is  declared  to  be  the  best. 
*3 


it  is  declared  to  be  the  best. 

\ 


Event 

{T  < 

X2 

Tv 

j  *2} 

(T  <T  , 

3 

•  •  •  •  •  •  • 

j  f»3} 

•  •  •  • 

(T  < 
*1- 

V 

j  f  k} 

Correct  decision 


Description:  ir^.  is  declared  to  be  the  best. 
Event ;  {T  <T  for  all  j  *  2,3,... ,k}. 

1  j 


Calculation  of  probabilities  of  wrong  and  correct  decisions. 

Let  b  *  cn,  g^(*)  the  probability  density  function  of  a  Gamma  random  variable 
with  scale  parameter  equal  to  unity  and  shape  parameter  equal  to  b  and  G^(-)  the 
distribution  function  associated  with  g^CO.  Note  that  nT^  has  Gamma  dis¬ 

tribution  with  scale  parameter  equal  to  unity  and  shape  parameter  equal  to  b.  Le 

"rWl-92 . V 

*  **ra  a  a  (committing  wrong  decision  of  type  r) 

9l'V’-',9k 

-  Pr  (T  <  T  for  all  j  + r) 

9l’e2**-*»9k  ir  ij 

*  P\.e2 . 0k<nTlr/e(r)  '  <nTij/9a))(9(J)/9(r))’  J,lr) 

■  .  0V2 . \  <nIi/9a) >  <9<r>'9<j>>*  f°r  a11  J 


since  T.  is  independent  of  T.  ’s,  j  i*  r. 
r  j 


WCWCWVISW 


^c«;v\w 


k 

!  n  (l- 
Jo  1-1 


Gb^0(r)/0(j))x^  8bW  dx* 


since  T^’s  are  independent, 
for  r  -  2,3,.. . ,k. 


P^ (®^»® 2* *  *  * ,0k* 

»  Prfl  .  (R  takes  the  correct  decision) 

“ 1*^2*  *  *  * 

-  Prfl  .  (T  < T  for  all  j  -  2,3,..., k) 

01,02’ *  *  * ,0k  H 

-  lo  "2(1-S((0(l)/0(j))x))  8b(x)  dx* 


(2.1) 


(2.2) 


Penalty  functions 

We  consider  two  types  of  penalty  functions. 

Penalty  functions  of  type  1 

Let  0-  { (0^,e2, . . . ,9^) ;0^ > 0  for  all  i}. 

For  every  (e^,92> . . . ,6^)  in  0  and  j-2,3,...,k,  let 

(0^,92,.  ..  ,®  ^) 

-  Penalty  for  taking  wrong  decision  of  type  j  when  the 
configuration  of  the  parameters  is  (0^,@2 . ,0^) 


c2  l08(9(2)/9(l)> 


for  some  c2  >  0  which  is  measured  in  the  same  units  as  those  of 

c  and  c. . 
o  1 


1 


The  rationale  behind  these  penalty  functions  is  explained  as  follows.  A 
critical  point  in  the  analysis  of  selection  problems  is  the  ability  to  discrimi¬ 
nate  between  the  best  population  it  with  the  parameter  value  6...  and  the  second 
best  population  tt^  with  the  parameter  value  ®^2)*  we  accePt  second  best 

population  as  the  best,  the  penalty  is  C2  logCB^/®^) •  The  penalty  for  other 
types  of  wrong  decisions  should  be  at  least  the  one  as  the  above,  and  for  mathematical 
expediency,  we  take  this  penalty  to  be  exactly  the  same  as  above.  The  loss 
function  now  works  out  to  be 


I>(® 2. » ® 2  *  •  •  •  > ®j^» 

*  Cost  of  sampling  +  expected  penalty 
k 

-  C(n)+c2  ^y.^(log(0 (2)^®(i)^j . 6k;  n) 

for  all  (01»®2 . ®k)  it*  ®  and  n  in  where  N  is  the  set  a21 

natural  numbers . 


(Ofosu  (1972)  has  taken  Wj  (®1»®2»  *  *  *  »\)  to  be  e9ual  to  c2  lo8^0(j)^0(l)^  *  j“2,3,.. 
This  is  more  general  than  the  one  we  considered  above.  Mathematically,  with  this 
choice  of  penalty  functions,  the  loss  function  seems  to  be  intractable  to  carry 
out  optimization.) 

The  loss  function  introduced  above  is  defined  over  the  cartesian  product  space 


0  x N.  In  order  to  find  the  minimax  sample  size  n,  we  maximize  L(e^,62,. . . ,6^;  n) 
over  all  (9i»92** *  * ,0k^  *n  ®^or  evet7  fixed  n  in  N.  Then  the  minimax  sample  size 
is  obtained  by  minimizing 

max  LCQjn) 

9e0  “ 


over  all  n  in  N.  We  will  take  up  this  work  in  Section  3 
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% 

V"v 


Penalty  functions  of  type  2^ 

Let  0  <  <S  <  1  and  a>0  be  prescribed.  The  constant  a  is  measured  in  the  same 


units  as  those  of  c  ,  c.  and  c„.  We  consider  the  following  penalty  functions. 

o  1  L 


For  (61,82,...,0jt)  in  0  and  j  »2,3,...,k,  let 


wj(e1,e2 . ek)  -  a  if  0(1)/0(j)  <6, 


0  if  6(1)/6(J)  >  6. 


The  rationale  behind  these  penalty  functions  is  as  follows.  Let  e^»92,****9k 
be  a  configuration  of  the  parameters  of  the  populations  iri»ir2****  »irjc  respectively. 


Suppose  the  statistical  procedure  R  declares  ir  to  be  the  best  population  for  some 


j 


j  a2,3,...,k.  The  parameter  values  associated  with  the  best  population  u  and  the 

11 


population  ir^  are  8^  and  8^  respectively.  If  the  ratio  C^.1)  of 


8^  and  8^j  is  close  to  unity,  we  would  not  like  to  be  penalized  for  taking  the 


wrong  decision  of  accepting  n  as  the  best  population.  On  the  other  hand,  if  the 


ratio  is  we  certainly  wish  to  be  penalized  for  accepting  ir^  as 

the  best.  A  line  has  to  be  drawn  somewhere  between  the  statements  that  the  ratio 


8 (1) / 8 ( j >  being  close  to  unity  and  that  it  is  being  small.  The  number  <$  distin¬ 


guishes  these  two  statements  and  the  choice  of  6  is  subjective.  We  assume  that 


the  penalty  for  accepting  the  population  tt  (j  »2,3,...,k)  as  the  best  when  the 

i 


ratio  9(i)/0(j)  *s  afflall»  i.e.,  <5  ,  is  the  same  constant  a.  The  loss  function 
then  works  out  to  be 


cv 

L^l’^2*  *  *  *  *®k*  (0^ » ®2*  *  *  *  * S^Pj  (8]_* ® 2 *  *  *  *  *  ^k*n^ 


for  (e1,02,...,0k)  in  0  and  n  in  N. 


For  a  given  configuration  C0]_»02’ *  *  *  *9k^  in  0,  one  could  explicitly  calculate  the 


loss.  In  order  to  find  the  minimax  sample  size,  we  need  to  maximize  the  loss 
function  over  the  parameter  space  Q  for  every  fixed  n  in  N.  We  will  take  up  this 
problem  in  the  next  section. 

3.  Minimax  Sample  Size 

In  this  section,  we  maximize  the  loss  function  L(*:  n)  over  the  parameter 
space  0  for  every  fixed  n  in  N  under  each  of  the  penalty  functions  of  types  1 
and  2.  First,  we  take  up  the  case  of  penalty  functions  of  type  1.  The  following 
subset  of  0  plays  a  crucial  role  in  the  maximization  problem  of  the  loss  function. 

0o  -  {(01,02,....ek)  e  0;  0X  -  r'V  02=  e3“  -  •  •  “  ek 

for  some  r  >_  1 } . 

(In  the  literature,  0q  is  the  so-called  set  of  least  favourable  choices.) 

The  following  chain  of  results  helps  to  evaluate  the  maximum  of  the  loss  func¬ 
tion  over  0  explicitly. 

Lemma  3.1  For  each  (0^»02*  *  *  *  * 9k^  e  ®o  w^t*1  91  “  r  ^92* 

P1C01,9?,...,0k;  n)  -  1°° (l-Gb(x/r))k  ^(xjdx, 

and 

PjC91,02,...,0k;  ti)  =|°° a-Gb(xr))Cl-Gb(x))k~2gb(x)dx 
for  every  j  *2,3,...,k. 

Further, 

max  p, (0jn)  ■  max 

0e0  1  r>l 

=*  o  — 

Proof.  By  direct  verification.  (Use  (2.1)  and  (2.2).) 


(l-Gb(x/r))k  1gb(x)dx. 
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» 


w  - 

P: 

$ 


Lemma  3.2  Let  (.9  ^ » 9  £ » •  •  •  *6^)  £0*  Write  8^  =  r  (Obviously,  r>_l.) 

•jf  tjf  ff 

Let  9,  *  9  ,,  v  and  9*0=...*  0,  *9,...  Then  (9.  ,9« , . , .  ,9.  )  e  0  and 
1  (.1)  i  3  k  \l)  1  l  Ic  o 


0  2*  *  *  *  *  ^k’  — ?1  ’  0  2  *  *  *  *  *  ^k’ 


Proof.  From  (2.2) 


P^( .9 »® 2 *  *  *  *  * = 


'o  j;2a-GbC(ea)/9(j))T),gb(lt)d!t 


r* 


(l-Gb(x/r))k  1gb(x)dx, 


since  1/r  -  9 (1)/0 (2)  i0(1)/0(3)  £•••  -0(l)/e(k)# 

is  is  is 

It  is  obvious  that  (0.  ,0„ , . . .  ,0,  )  £  0  .  By  Lemma  3.1, 

1  L  k  o 


*  *  *  r° 

p^(9^j02» • • • n)  *  J  (1  — 


Gfe(x/r))k  1gb(x)dx. 


This  completes  the  proof. 

The  following  theorem  simplifies  the  problem  of  maximization  mentioned  earlier. 
Theorem  3.3  For  every  n  in  N, 


,w 

1>V 


.1 

I 

'4 

d 


N.i 

ty 

1 


max  LC§j  n)  *  max  L  (Jj  n) 

0E0  9e0 


Proof  Let 


and 


C(n)+c,  max  Clog  r)  (1  -  J  (1-G  (x/r))k  1g  (x)dx) 
z  r>l  >0 

roo 

C(n)+c2(k-l)  max  (log  r)  j  (1  -  Gfa(xr) ) (1  -  Gb(x))k  2gb(x)dx. 


max  l  (log -r-^-)p  (0  ,9  ,...,9,  ;n) 
0£0  j=2  y(l)  J  1  1  K 


k  ^  (2 )  is  is  is 

max  l  (log  -^-)p  (0  ,0  ,...  ,9  ;n). 

0*£0  j-2  0...  J 

=  (1) 
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To  prove  the  first  part  of  the  theorem,  it  suffices  to  show  that  p*q.  Since 
0Qe0,  q^p.  We  now  prove  that  p£q.  For  this,  it  suffices  to  show  the  following 

if  if  if 

Given  d.^,02****  »®k)  e  t*1®1®  exists  (0^,02»  ••  •  »9k)  e  ®0  suc^  c^at 
k  0 


y.  dog  T^Ph  C01*09,...,0lr ;n) 
J-2  0d)  J  1  2  k 


k 


£  l  dOg  ^  )  P  .  (6  -I  »  S  r\  »  •  •  •  »  O)  , 

-1=9  Q  J  L  4  K- 


*  * 


j-2 


(1) 


*  * 


Let  Sj-I(1),  VV’" 'VW 


*  * 

Obviously,  “r  9(l)  ^or  some  an^ 


if  if  if 

C® €  0  .  Also , 

1  Z  K  O 


I  iogCo  )p.  C0 ,  »9«  »  •  •  •  ,9,  ;  n) 

j-2  (1)  J  1  ^  * 


=  ^  (log  r)  p  (0  ,09,...,0  ;n) 

j-2  j  1  2 


dog  r)  (i-p1(91,02 . ®k;n) 


*  * 


1  (logr)  (1  "P]/6].’  92».*»,9k;n), 


by  Lemma  3.2 
k 


1  (logr)  l  p,(e*,e*,...,e";n) 

j«*2  J  L  *  K 


y  log(-^-)p  (0  ,0  , ...  ,0.  ;  n) 
W  e(1)  J 


*  * 


This  shows  that  paq  and  hence  max  L(^;n)  *  max  L(,0;n).  Lemma  3.1  provides  the 

0e0  0e0 


exact  expressions  for  max  L(9;n) . 

0E0 

=  O 

In  view  of  the  above  theorem,  the  critical  function  in  the  determination  of 


the  minimax  sample  size  is  the  following  one: 


f  Cr)  -  (log  r)  (1  - 


k-1 


(l-Gb(x/r))  gfe(x)dx,  r^l. 


where  b  ■  cn. 
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We  undertake  the  study  of  the  above  function  in  Section  4.  Now,  we  work  out 
the  niay-ftnum  of  the  lass  function  L^yn)  under  the  penalty  functions  of  type  2. 

Let  0  <  6  <  1  and  a  >  0  be  fixed.  We  want  to  maximize  L^n)  over  all  £e  0 
for  each  sample  size  n  under  the  penalty  functions  of  type  2.  See  Section  2. 

For  this  purpose,  we  introduce  the  following  subset  of  0. 

0^  “  *"  *  ®k^  ’ 

The  following  theorem  works  out  explicitly  the  desired  maximum  of  the  loss 
function. 

Theorem  3.4  The  following,  statements  are  true. 

(1)  If  (01,e2,...,0k)  e0lf  then 

p^  (®2.»  ®2 *  *  *  *  *  =  j^0 ^  ”  ®{)  ^  g^(x)dx. 

C2)  If  (QL,e2,...,9kl  e0lt  then 

1>(6^»92»»«*  *®k*n^ 


n. 

C(n)  +a  l  p.(9.»9,,...,6  ;i 
j-2  j  1  1 


■  CCn)  +  a(l  “Pi(®i»®2**#**®k’n^' 


(.9^,92*  •••  ,9k)c0 


*JC®i»®2***  *  *®k*n^ 


L(9^,92,...  »©jtJn) 


C®1*®2*  *  *  * *®k^e®l 

C(n)+aCl-|" Cl  -Gb(6x))k"’1gb(x)dx). 


(Recall  that  b*cn.) 
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Proof. 

(1)  If  ce1,e2>...,ek)  e01,  then  |  -02  -03  - ... -8k.  So,  0(1)-01, 

9(2)  “92,***’9(k)  *8k  accordlnS  t0  our  convention»  Therefore,  from  (2.2), 

W92 . 6k!n)  '  fjS2Ci-Gb<;9^x))8b<!c)d!t 

o  ^ 


r°° 
.  0 


(1  -  Gb(6x))k-1gb(x)dx. 


(2)  If  (-01»62 . 0k)e0lt  then  (01,02,...,0k)  -  a  for  all  j-2,3,.. 

Consequently, 


•  )k« 


L(01,02,...,0k;n)  -C(n)  +^2Wj(01,02 . 9k)Pj (ei*e2* *  *  * ,0k;n) 

k 

“  C(n)  +  a  £  P  (9.  ,0_ , . .  •  ,8,  ;n) 
j-2  3 

-C(n)  +a(l-P1(91,92,...,0k;n)) 

-  C(n) +a(l-|°° (l-GhCS*))11  1gb(x)dx), 
by  (1)  above. 


(3)  Since  0^0, 


p  *  max  L(0. ,9_, . . . ,9.  ;n) 

C0lf02,...,ek)e0 

>_  max  L(01 ,9-,...  ,9.  ;n)  «q,  say. 

(01,02,...,0k)e01 

r°  k-i 

From  (2)  above,  q«CCn)+a(l-J  Cl-Gb(6x))  g^(x)dx).  To  complete  the  proof,  it 
is  enough  to  show  that  p<_q.  For  this,  we  show  that  for  any  (03,02,...  ,9k)  in  0, 

LC01,92,...,8k;n)  <_q. 
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Let  (91,92,...,0k)  e  0.  Let  0^10^1  •••l9(k)  be  the  ordered  arrangement  of 
0-.  ,0.,. . .  ,0,  In  Increasing  order  of  magnitude. 

■L  4 

Case  (i).  ^  >  5  for  all  j  -2,3 . k.  Then 

-  9«) 


L(®1*02**** *®k*n)  “  C(tl)+01.<J»  Obviously. 

®(1)  * 

Case  (ii) .  There  exists  j  in  {2,3,...  ,k}  such  that  Let  j  be  the  smallest 


integer  in  {2,3,...,k>  such  that 


Ii)_< 

'(j*) 


(j) 

6.  This  means  that 


and 


Then 


(i) 


6  for  all  i- l,2,...,j* -1 


-^•<5  for  all  i  -  j*,  j*  +  l,...,k. 
6(i)  " 


L(01,02,...,0k;n)  -  C(n)+a  J  Pj^.Oj . 0fc; 

k 

<_  C(n)  +a  £  P  (0_,0_,...  ,0,  jn) 
j-2J 


n) 


<  C(n)  +a(l -P1(01,e2,...,0k;n)). 


Let  us  calculate  P^(0^,02>... ,0k;n).  Note  that 


and 


Therefore, 


*^•<1  for  J  -2,3 . j*-l 

6<J) 

<  t  for  j  -  j*,J«+l . k. 

9CJ)  ~ 
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■  W  x  £  x  for  j  *  2,3,...  ,j*-l 

9(j) 

and 

6 

x <  Sx  for  j 
9(j)  " 

for  all  positive  x. 

Consequently, 


Pl(0l,82’***,9k;n) 


0 

a  Cl-G.  <  -J&  rt)t(x)4« 
J-2  b  eU)  ^ 


x))gb(x)dx 


>  j" Cl-Gb(x))J*"2 


(1  -Gb(6x))k“;J*+1gb(x)dx. 


From  the  above  inequality,  we  obtain 
L(0^,02»»»» » 

<  C(n)  +a(l-P1(01,02,...,0k;n)) 

£  C(n)  +aU-J  (l-Gb(x))J*"2(l-Gb(6x))k';1*+1gb(x)dx). 

Now,  we  show  that 

j" (1  -  (x) ) j  *”2 ( 1  -  Gb ( 6x) )k”j  *+lgb (x) dx 

>  j" (1 -Gb(6x))k”1gb(x)dx. 

Since  0  <  <5  <  1,  G^iSx)  <Gb(x)  for  all  positive  x  from  which  it  follows  that 
Cl -Gb(x))j*”2  <  (1  -GhCSx))-5*”2  for  all  positive  x.  This  inequaltiy  yields 


'UT  vnr’tnr  cn 
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L(01,02,...,0kJn) 

r 

<_  C(n)  +a(l  -  I 

Jo 

<  CCn)  +  aCl-|  (1  - 

J  0 

This  completes  the  proof. 

In  the  case  of  penalty  functions  of  type  2,  the  maximum  of  the  loss  function 
L(J3,;n)  over  0  works  out  quite  explicitly.  In  order  to  find  the  minimax  sample 
size,  we  need  to  minimize  this  maximum  with  respect  to  n.  We  take  up  this  problem 
in  Section  4. 


GbC5x))k”1gb(x)dx)  -q. 


(1  -  Gb(x))j*"2Cl  -  Gb(«Sx))k“j*+1gbCx)dx 


4.  Numerical  Illustrations 

The  selection  of  the  best  Gamma  population  has  considerable  bearing  in  the 
selection  of  the  population  which  is  least  variable.  Let  . . .  ,irk  be  k  normal 

populations  with  i-th  population  having  mean  and  variance  0i  >  0.  The  best  popu¬ 
lation  among  these  is  the  one  with  the  smallest  9^.  If  X^j ,  j  *1,2 . n  is  a 

ramdom  sample  of  size  n  from  it.  ,  then  we  could  use  the  statistics  T.  *  £  CX  ”x.«)  » 

j-1  i;J 

i«l,2,...,k  to  discriminate  the  populations  »^k»  T^/2B1  has  Gamma 

distribution  with  shape  parameter  equal  to  (n-l}/2  and  scale  parameter  equal  to 
unity.  In  this  case,  in  order  to  find  the  minimax  sample  size,  one  has  to 
maximize  ' 

LCrjnl  -  C(n)  +C2Clog  r)  Cl¬ 
over  r  ^1  for  every  fixed  n»  -,3,...,  and  then  this  maximum  is  to  be  minimized  over 
the  set  {2,3,4,...},  where  b  *  (n-l)/2.  In  this  section,  we  address  ourselves  to 
the  problem  of  determining  minimax  sample  sizes  when  b  *  n/2. 

The  critical  function  in  the  maximization  of  L(r;n)  over  r  >_1  is  the  function 

*njk(r)  •  log  rCl  -  j^Cl -Gb(x/r})k“1gb(x)dx),  r>_l. 


Cl  -GbCx/r))k_1gbCx)dx) 
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where  b  -  n/2.  Determination  of  the  maximum  of  this  function  analytically  does 


not  seem  to  be  feasible.  It  is  obvious  that  f  .  (1)  -  0  and  one  can  prove  that 

lim  f  (r)  "0.  An  extensive  tabulation  of  this  function  over  the  domain  of 
_  n«K 

jMoo  9 


definition  [1,*)  of  r  for  various  values  of  n  and  k  suggests  that  this  function 


is  unimodal.  These  calculations  yielded  the  following  tables.  (Let  r  ,  denote 

the  number  at  which  f  ,  (r  ,  )  is  maximum.) 

n,k  n,k 


TABLE  1 


No.  of  populations 
k 

2 

3 

4 

5 

6 

sample  size 
n 

rn,k 

rn,k 

rn,k 

rn,k 

Tn,k 

2 

4.06 

5.21 

6.26 

7.24 

8.18 

3 

3.09 

3.77 

4.34 

4.85 

5.32 

4 

2.69 

3.20 

3.60 

3.96 

4.27 

5 

2.47 

2.88 

3.21 

3.48 

3.72 

6 

2.32 

2.68 

2.96 

3.19 

3.39 

7 

2.22 

2.54 

2.79 

2.99 

3.16 

8 

2.15 

2.44 

2.66 

2.84 

3.00 

9 

2.09 

2.36 

2.57 

2.73 

2.87 

10 

2.04 

2.30 

2.49 

2.65 

2.77 

11 

2.00 

2.25 

2.43 

2.57 

2.69 

12 

1.97 

2.21 

2.38 

2.52 

2.63 

13 

1.95 

2.18 

2.34 

2.46 

2.57 

14 

1.92 

2.14 

2.30 

2.42 

2.52 

15 

1.90 

2.12 

2.27 

2.39 

2.48 

16 

1.89 

2.10 

2.24 

2.35 

2.45 

17 

1.87 

2.07 

2.22 

2.33 

2.41 

18 

1,86 

2.06 

2.20 

2.30 

2.39 

19 

1.85 

2.04 

2.18 

2.28 

2.36 

20 

1.84 

2.03 

2.16 

2.26 

2.34 

The  following  pattern  emerges  among  r  ,  *s  from  the  above  table. 

n,«t 


(i)  r  ,  <  r  ,  if  m  >  n. 

m,k—  n,k  — 

(ii)  r  ,<r  if  k<s, 

ntk—  n,s  — 


Using  the  property  (i)  mentioned  above,  we  promulgate  the  following  strategy 
to  determine  minimax  sample  sizes. 

Strategy,  k  is  fixed.  Find  r^  the  point  at  which  f^  k(r)  is  maximum.  Evaluate 
L(r;n)  for  r-1.0,  1.1,  1.2,, ...r^  ^  and  n  ■  2,3,...  ,20.  These  values  are  tabulated 


in  a  two-way  grid,rows  corresponding  to  r  and  columns  to  n.  For  each  column  in 
the  two-way  grid,  locate  the  maximum.  Under  each  column  representing  a  particular 
sample  size  n,  we  know  that  the  maximum  of  L(r;n)  occurs  at  some  value  of  r  between 
1.0  and  r.  .  .  See  property  (i)  above.  Then  the  column  for  which  this  maximum  is 

2  f  K. 

minimum  is  located.  The  corresponding  sample  size  is  the  required  minimax  sample 
size. 

If  the  column  maximums  are  increasing  right  upto  r  *20,  then  the  minimax  sample 
size  is  ^20.  In  such  a  case,  we  need  to  include  some  more  sample  sizes  n* 21,22,... 
etc.  For  largevalues  of  n,  evaluation  of  L(r,n)  may  be  beyond  the  reach  of  many 
computers.  Later,  in  this  section,  we  give  asymptotic  formulas  for  L(r;n). 

The  above  method  of  determining  minimax  sample  sizes  is  a  modification  of 
a  procedure  suggested  by  Ofosu  (1972).  This  modification  works  faster  than  the 
methods  outlined  by  Ofosu  (1972). 

By  way  of  illustration  of  the  working  of  the  above  strategy,  we  include 
the  following  examples. 

1.  No.  of  populations,  k“2,3,4,5,6. 

2.  Cq  *  0,  c^  -  1,  d  -  1  so  that  C(n)  -  kn. 

3.  c2  - 50(25)  500. 

The  loss  function  then  becomes 

•00 

L(r;n)  *kn+c2(log  r)  (.1  -  (1  -Gb(x/r))  "^(xidx), 

r>_l,  n  ■  1,2, 3, 4,. . . ,  where  b»n/2. 

The  minimax  sample  sizes  are  given  in  the  following  Table. 
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Table  2 

Minimax  sample  sizes  (n)  and  minimax  loss  L(n)  under  penalty  functions  of  type  1. 


_ _ id 

_ l _ 

< 

4 

5  1 

6 

C2 

n 

LCn) 

n 

,  L(n) 

n 

'LCn) 

n 

•  LCn) 

n 

L(n) 

50 

3 

27.79 

3 

42.60 

3 

53.79 

3 

63.07 

3 

71.18 

75 

4 

38.55 

4 

52.25 

4 

72.87 

4 

84.90 

4 

95.63 

100 

4 

48.74 

5 

73.23 

5 

90.99 

5 

105.53 

5 

118.15 

125 

5 

53.76 

5 

87.79 

6 

108.69 

6 

125.63 

5 

140.19 

150 

5 

68.53 

6 

101.97 

6 

125.63 

6 

144.75 

6 

161.15 

175 

5 

78.29 

6 

115.97 

7 

142.57 

7 

163.80 

6 

182.01 

200 

6 

87.95 

7 

129.79 

7 

158.93 

7 

182.20 

7 

202.02 

225 

6 

97.44 

7 

143.39 

7 

175.30 

7 

200.60 

7 

222.02 

250 

6 

106.94 

7 

156.99 

8 

191.33 

8 

218.47 

8 

241.58 

275 

6 

116.43 

8 

170.43 

8 

207.26 

8 

236.32 

8 

260.94 

300 

7 

125.90 

8 

183.74 

8 

223.20 

9 

254.05 

9 

280.17 

325 

7 

135.23 

8 

197.05 

9 

238.82 

9 

271.47 

9 

299.02 

350 

7 

144.55 

9 

210.36 

9 

254.42 

9 

288.90 

9 

317.86 

375 

7 

153.88 

9 

223.45 

9 

270.02 

10 

306.11 

10 

336.56 

400 

7 

163.20 

9 

236.55 

10 

285.46 

10 

323.18 

10 

355.00 

425 

7 

172.53 

9 

245.65 

10 

300.80 

10 

340.25 

10 

373.43 

450 

8 

181.76 

9 

262.74 

10 

316.14 

10 

357.33 

11 

391.85 

475 

8 

190.96 

10 

275.77 

11 

331.36 

11 

374.18 

11 

409.95 

500 

8 

200.17 

10 

288.71 

11 

346.48 

11 

390.97 

11 

428.05 

Nov,  we  consider  Che  case  of  penalty  functions  of  type  2.  From  Theorem  3.4, 


max  L(0_:n) 

S.60 


COO  +a(l  -  j" (.1  -G^CSxJJ^^CxJdx, 


for  n  ■  1,2,3,...,  0<6  <1  and  b  -  cn.  In  order  to  find  the  minimax  sample  size, 
we  need  to  minimize  the  above  function  over  n-1,2,3,...  for  a  given  k,<5,  c  and 
C(n)  *  kn.  As  an  illustration,  we  cake  c  ■  --  and  a  » .100(100)600  and  determine  the 


L9 


i 

Table  3 


i 


I 

6=0.50 


I 

i 


:  6=0.60 


,  6=0.70 

» 


i 

! 

6=0.80 

i 


Minimax  sample  sizes  (n)  and  minimax  loss  function  L(n)  under  penalty  functions 


of  type  2. 


Table  3  (continued) 

Minimax  sample  sizes  (n)  and  minimax  loss  function  L(n)  under  penalty  functions 
of  type  2. 


\k 
a  \ 

n 

L(n) 

V 

n 

\ 

L(n) 

4 

n 

L(n) 

C 

n 

> 

L(n) 

100 

2 

51.56 

2 

70.48 

2 

81.17 

2 

88.46 

200 

2 

99.13 

2 

134.96 

2 

154.33 

2 

166.91 

300 

3 

146.30 

2 

199.44 

2 

227.50 

2 

245.37 

400 

7 

192.21 

3 

263.69 

3 

301.68 

3 

326.10 

500 

7 

236.76 

1  4 

327.34 

3 

374.10 

2 

402.27 

600 

9 

■  281.05 

i  7 

389.15 

.  3 

446.52 

.  2 

480.73 

5.  Asymptotic  Minimax  Value  of  a 

In  this  section,  we  approximate  the  following  integral 


C 


Cl-GnCf))k"1«n(*)dx 


for  large  values  of  n  by  an  integral  involving  the  standard  normal  distribution, 

where  g  (•)  is  the  probability  density  function  of  a  Gamma  distribution  with 
n 

parameters  n  and  1,  i.e.. 


*»<*>  ■  rfer  0<x<-t 


and  G  (•)  is  the  distribution  function  associated  with  gn (* ) •  This  integral  appears 
a  n 

in  the  loss  functions  considered  under  penalty  functions  of  types  1  and  2.  This 
integral  can  be  realized  probabilistically  as  follows. 

Let  Y. ,  Y, , . . .  ,Y,  be  k  independent  identically  distributed  random  variables 
with  common  probability  density  function  gnC»)  given  above.  Then  for  any  r  >.1, 

Pk(r;n) -Pr(Y1<  rY  ,  j-2,3,...,k) 


r « 

h  j-; 


PrCrYj  >  x)  gn(x)dx 
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We  observe  chat  2Y^''^x2n  for  every  j*l,2,...,k.  We  want  to  apply  Wilson- 
Hilferty's  transformation  to  each  2Y  ^ .  See  Kendall  and  Stuart  (1977,  p.  398-399). 


Let 


^2n  ’  u  9n' 

(— )1/2 
v'9ni 


for  j  *  1,2,... ,k. 

Zi,Z2,. •• ,2^  are  independent  random  variables  and  for  large  n,  each  Z^  has  standard 
normal  distribution.  We  also  note  that 


2Yl  <  2zYj 


if  and  only  if 


TEpTT 

'k9nt 


<-V« 


r1^3  [  ( — 1)  1^3  ^ _ 1_)  ] 

u2n  '  u  9n'J 

cX?75 

(1  -  ~)Cr1/3  -  1) 

frT772 

v9n' 


if  and  only  if 


for  any  j  ■  2,3,... ,k,  where 


7  <  r-ll37  -L  a 

4.  ^  r  l.  +  a 
1  j  n,r 


Cl-£Xrl/3-l> 


a 


n,r 


Tm75" 

v9n; 


Now,  we  are  in  a  position  to  approximate  the  integral  in  focus  as  follows. 
For  large  n. 


rvi 


^Cl-Gtt(|))lt'1giiU)<lx  -  Pk(r;n) 

-Pr(^<rYjf  j  -2,3,. ...k) 

•  Pr  (2YX  <  2rY^  ,  j  « 2,3,... ,k) 

1/3 

■Pr(Z.  <rA/JZ. +a  ,  j-2,3,...,k) 
l  J  n,r 

-  P  (.1  -*( - - ^Ypi  ®*PC-  ^-)dx» 

J—  r  (2it) a/  z 

where  $(•)  is  the  distribution  function  of  the  standard  normal  probability 

model. 

The  loss  function  under  penalty  functions  of  type  1  for  alrge  n,  becomes 

C  X  *  a 

L(r;n) -C(n)+c2(logr)  (1-1  (l-*( - 

*  •*  r 


"17?  «PC-  — )d*) 


for  r>_l  and  n  in  N,  where  b  »  cn. 


To  find  the  minimax  sample  size  n,  the  above  approximation  can  be  used  in 
the  strategy  outlined  in  Section  4. 

The  maximum  of  the  loss  function  over  the  entire  parameter  space  under 
penalty  functions  of  type  2,  for  large  n,  becomes 


C(n)  +a(l 


-j>  % 


bd.«5_k-l  1  2  .  . 

YT7l-))  - 172  e  dx) 

Cy)  (2iT)i/Z 


for  n  in  N,  where  d  *  cn  and 


a  -  £>«£) 1/3  - 1) 

^9d; 


Similar  approximations  can  be  provided  for  the  problem  of  selecting  the 


best  normal  population  with  least  variance. 


f « 
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